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ABSTRACT 

We present 16-GHz Sunyaev-Zel'dovich observations using the Arcminute Microkelvin Im- 
ager (AMI) and subsequent Bayesian analysis of s ix galaxy clusters a t redshift {z ~ I) chosen 
from an X-ray and Infrared selected sample from lCulverhouse et alJ (12010 ). In the subsequent 
analysis we use two cluster models, an isothermal yS-model and a Dark Matter GNFW (DM- 
GNFW) model in order to derive a formal detection probability and the cluster parameters. We 
detect two clusters (CLJ1415+3612 & XMJ0830-I-5241) and measure their total masses out to 
a radius of 200 x the critical density at the respective cluster's redshift. For CLJ1415H-3612, 
we find Mt,2oo = 7.27+}-^^ x IO'^^Mq (J3 model) and Mj.aoo = 10.402f3g x IO'^^Mq (DM- 
GNFW model) and for XMJ0830-I-5241, we find Mt,20() = 3.56![ J'^ x'iO^'^Mq, (yj-model) 
and Mt,20o = 4.66+| 4j x 1O''*M0 (DM-GNFW model), which agree with each other for each 
cluster We also present maps before and after source subtraction of the entire sample and 
provide ID and 2D posterior marginalised probability distributions for each fitted cluster pro- 
file parameter of the detected clusters. Using simulations which take into account the mea- 
sured source environment from the AMI Large Array (LA), source confusi on noise, CMB pri- 
mordi als, instrument noise, we estimate from low-radius X-ray data from lCulverhouse et al.l 
(1201 Ol) the detectability of each cluster in the sample and compare it with the result from the 
Small Array (SA) data. Furthermore, we discuss the validity of the assumptions of isother- 
mality and constant gas mass fraction. We comment on the bias that these small-radius 
estimates introduce to large-radius SZ predictions. In addition, we follow-up the two de- 
tections with deep, single-pointed LA observations. We find a 3 - cr tentative decrement 
toward CLJ1415H-3612 at high-resolution and a 5 - cr high-resolution decrement towards 
XMJ0830+524L 

Key words: cosmology: observations - cosmic microwave background - galaxies: clusters 
- Sunyaev-Zel'dovich X-ray- galaxies: clusters: individual (CLJ1415-H3612, ISCS1438-H34, 
RXJ0910-H5422, SPJ1638+4039, XMJ0830+5241, XMJ0849-H4452) 
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1 INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zeldovich 1970, 
1972) is the inverse-Compton scattering of the CMB radiation by 
hot, ionised gas in the gravitational potential well of a cluster of 
gala xies; for reviews on the SZ effect see e.g. iBirkinshawl h999l) 
and ICarlstrom, Holder, & Rees3 ( l2002h . The effect is useful in a 
number of ways for the study of galaxy clusters; here we are con- 
cerned with two in particular. First, because the SZ effect arises 
from a scattering process, a cluster at one redshift will produce the 
same observed SZ surface brightness as an identical cluster at any 
other redshift. On the other hand, the integrated SZ flux does de- 
pend on the angular diameter distance and therefore the redshift 
z, but the dependence is small at z g 0.5. Hence, the usual sensi- 
tivity issue of high-redshift observing is avoided, which is particu- 
larly useful for this study. Second, since the SZ surface brightness 
is proportional to the line-of-sight integral of pressure through the 
cluster, the SZ signal is less sensitive to concentration than the X- 
ray Bremsstrahlung signal; one corollary of this is that the ratio 
SZ-sensitivity / X-ray-sensitivity increases with distance from the 
cluster centre so that with SZ one can probe out to, say, the virial ra- 
dius, provided the telescope is sensitive to sufficiently large angular 
scales. SZ decrements are faint, however, and can be contaminated 
or obliterated by sources of radio emission. Many SZ studies are 
currently and routinely carried out at different redsh ift ranges by 
the Atacama Cosmology Telesco pe (ACT) (se e e.g. Sehgal e t al. 
2012), the Planck satellite (see e.g. lPlanck and A MI Collabor ations! 
2012h . the South Pole Telescope (SPT) (see e.g. lstalder et al.l20i; 



Table 1. AMI technical summary. 



and t he Sunyaev-Zel'dovich Array (SZA, see e.g. Muchovei et al.l 
l2007h . The Arcminute Microkelvin Imager (AMI) has demon- 
strated its capability of stu dying X-ray se lected clusters of galax- 



ies at moderate z, see e.g. lAMI Consortium: J gwar t et al. (2011) 



IaMI Consortium: Rodriguez-Gonzalvez et all ( 1201 II) , however al- 
ways at relatively low redshift. Here we investigate the feasibility 
of observing and detecting clusters at higher (z ~ I) redshift as well 
as constraining their masses. 



1.1 Sample selection 

ICulverhouse et al.1 ( 120101) selected a sample of 1 1 clusters a.t(z ~ 1) 
from infrared and X-ray studies, and observed these with the SZA 
to search for SZ signal. They found SZ decrements towards 3 of the 
clusters; using X-ray based values of ^2500, they al so estimated 72500 
and_Mj^j_250o (from the Y-Mg^s scaling relation of lBonamente et aP 
I2OO8I) . Here, we present AMI observations of the subsample of 
6 (of the original 11) clusters which are at 5 > 20°, along with 
Bayesian inference where possible. Our initial aims were to: 

(i) follow up the SZA observations with an instrument with dif- 
ferent systematics; 

(ii) consider the effects of difference in analysis procedures used 
by ICuIverhouse et al] ( I2OIQ) (from X-ray measurements near the 
cluster centre, and assuming hydrostatic equilibrium; these as- 
sumptions are often made in the literature) and by ourselves (tem- 
perature from SZ observation plus M-T relation from virial theo- 
rem, as well as Bayesian and probability search). 

The names, coordinates and redshifts of the sample are shown in 
Tabled 

In section |2] we briefly review AMI and then describe the 
observational programme in section [3] The data reduction and 
Bayesian analysis are explained in section|4] We have conducted re- 
alistic AMI simulations using X-ray measured parameters to model 





SA 


LA 


Antemia diameter 


3.7 m 


12.8 m 


Number of antennas 


10 


8 


Baseline lengths (current) 


5 -20m 


18- 110m 


Primary beam at 15.7 GHz 


20'. 1 


5' .5 


Synthesized beam 


a 3' 


« 30" 


Flux sensitivity 


30mJy s-'/2 


3 mjy s-"^ 


Observing frequency 


13.9- 18.2GHz 


13.9- 18.2GHz 


Bandwidth 


4.3 GHz 


4.3 GHz 


Number of channels 


6 


6 


Channel bandwidth 


0.72 GHz 


0.72 GHz 



Table 3. Noise levels of each tai'get on the central regions of the continuum 
maps and integration times on both SA and the corresponding LA raster 
observations. 



Cluster 
Name 


SA noise 

(pJy) 


LA noise 


tobs.SA 
(hours) 


tob.s,LA 
(hours) 


CLJ1415+3612 


95 


65 


36 


17 


XMJ0849+4452 


70 


80 


11 


22 


ISCS 1438+34 


90 


75 


71 


42 


RDJ09 10+5422 


100 


115 


34 


12 


SPJ1638+4039 


145 


85 


55 


18 


XMJ0830+5241 


70 


45 


44 


31 



the galaxy cluster and discuss their validity in section |5] present 
the SA analysis in section |6] and the LA follow-up of two clusters 
in section |7] Throughout, we assume Ho = 70 kms-'Mpc"' and 
a concordance ACDM cosmology with f2„,_o = 0.3, Q,\fl = 0.7, 
Hi = 0, fifc.o = 0.041, Wo = -I, w,, = and o-g = 0.8. 



2 THE ARCMINUTE MICROKELVIN IMAGER 

The Arcminute Microkelvin Imager (AMI) is a dual interferomet- 
ric array, near Cambridge. It consists of an 8-element Large Array 
(LA) with dishes of 12.8m and the Small Array (SA), a 10-dish 
array with an antenna size of 3.7m each. Both arrays observe at 
a central frequency of 15.75 GHz and their configuration is opti- 
mised for the specific scientific area they were conceived for. The 
main characteristics of A MI are summarised in Table[T]and are de- 
scribed in greater detail in lAMI Consortium: Zwart et al. I diooi) . 



3 AMI OBSERVATIONS 

All the clusters were observed using both arrays; a mosaicing strat- 
egy was used to cover the same area with the LA as the single 
pointed observations measured by the SA. Integration times and 
map noise levels for each target on both arrays are shown in Ta- 
ble[3] Some of the clusters had substantially more integration time 
(e.g. ISCS 1438+34) compared to others which reflects the amount 
of flagging necessary to remove interference and to assure a highly- 
filled circularly symmetric uv-coverage to sample all observable 
spatial scales and obtain a circular synthesised beam. All obser- 
vations were taken between May 2008 and March 2011 on the 
SA, while integration on the individual targets were all carried out 
within a few months to reduce the effect of source variability. Dur- 
ing each run on the SA we observed a close-by secondary cali- 
brator for 400 seconds each 6 minutes to maintain phase stability 
on the SA. The LA mosaicing runs were carried out within a few 
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Table 2. Cluster sample observed by AMI: redshifts, coordinates and the reference of the initial detection/redshift measurement 



Cluster name 


Redshift 


R.A. 


DEC. 


Detection 


CU1415+3612 


1.03 


14 15 11 


36 12 04 


X-rav. Maushan et al. (2006). Perlman et al. (2002) 


ISCS1438+34 


1.41 


14 38 09 


34 14 19 


IR. Stanford et al. (2005) 


RDJ0910+5422 


1.11 


09 10 45 


54 22 09 


X-rav. Stanford et al. (2002) 


SPJ1638+4039 


1.20 


16 38 52 


40 38 43 


IR, Muzzin et al. (2009^ 


XMJ0830+5241 


0.99 


8 30 26 


52 41 33 


X-ray, Lamer et al. (200^ 


XMJ0849+4452 


1.26 


08 48 59 


44 51 50 


X-rav, Rosati et al. (1999) 



days of each SA observation and each run had a secondary calibra- 
tor interleaved every 10 minutes for 2 minutes. Although the SA 
is our primary SZ array with baseline ranges geared to be more 
sensitive to SZ flux and larger scales, we followed-up two targets 
(CLJ1415+3612 and XMJ0830+5241) with deep LA single point- 
ings, as high redshift clusters are expected to have a smaller an- 
gular extent due to the angular diameter distance-redshift relation. 
SZ observations at high redshift could therefore benefit from the 
increased resolution of the LA. To properly assess this and get fur- 
ther detection confirmation, we carried out these deep follow-up 
observations on each target within a month in 2012. LA integra- 
tion time for the single pointings on CLJ1415-I-3612 was 8 hours 
in total and reached a sensitivity of 25 /jJy (50 yuJy on the shortest, 
< 1.5U basehnes). For XMJ0830+5241 a sensitivity of 20 yujy (45 
/jjy on the shortest, < 1.5k/l baselines) was reached after 13 hours. 
The pointed LA observations visit the phase calibrator every 10 
minutes for 100 seconds. 



4 DATA REDUCTION AND ANALYSIS 
4.1 Data reduction and mapping 

Individually, the data gathered from each run on both arrays are 
flagged for shadowing effects, slow fringe rates, pointing - and path 
compensator delay errors with our in-house reduction package re- 
duce. Absolute flux calibration is performed using daily observed 
calibrators, 3C48 and 3C286. The data are then Fourier transformed 
and fringe-rotated to the pointing centre. Further flaggings reject in- 
terference and discrepant baselines. Phase calibration is performed 
using the interleaved calibrators. Last, the amplitudes of the visi- 
bilities are weighted to take into account system temperature varia- 
tions, due to weather and airmass, before being outputted for every 
frequency channel as UVFITS. 

Our in-house software fuse is used to concatenate the UVFITS 
of each observation and then the data are mapped using the imag- 
ing package aip£] for both the single frequency channels and the 
full spectrum. An initial deconvolution to estimate the real map 
noise is done using the CLEAN-algorithm with a flux limit of 3 
times the noise calculated from the weights for each channel and 
pointing centre. The noise level of those initial maps are calcu- 
lated and serve as flux limit (3 x the calculated value) for a sec- 
ond deconvolution to produce the final maps. We use a box en- 
compassing the entire pri mary beam. For more details on the map- 
pin g technique, see e.g. AMI Con sortium: Franzen et al. 
and lAMI Consortium: Shimwell et alj ( l2012l) . 



http://www.aips.nrao.edu 



4.2 Source finding 

The LA maps are used to perform source finding and spectral in- 
dex fitting using our in-house software sourcefind. We give a brief 
su mmary of the technique here, for a mo re detailed description 
see lAMI Consortium: Franzen et al. 1 1 |201 ll) . All pixels on the map 
with a flux density greater than 0.6 x y x cr„, where cr„ is the 
noise map value for that pixel and y = 4 is the desired detection 
threshold, are identified as peaks. The flux densities and positions 
of the peaks are determined using a tabulated Gaussian sine de- 
gridding function to interpolate between the pixels and the peaks 
above a threshold of y x a„ are identified as sources. In addition 
the MPS routine jmfit is used to fit a two-dimensional Gaussian to 
each source to give the angular size and the integrated flux density 
for the source. These fitted values are compared to the point-source 
response function of the telescope to determine whether the source 
is extended on the LA map. 

Assuming a power-law relationship between flux density and 
frequency (5 oc v^") and given the measured flux of each individ- 
ual channel map, a spectra l index was calculat ed using an MCMC 
method called metromod ( iHobson et a"i]|2004l) — the prior on the 
spectral index has a Gaussian distribution with a mean of 0.5 and 
cr of 2.0, truncated at ±5.0. The fluxes and spectral indices are re- 
tained for use in our Bayesian analysis. 



4.3 Preparing SA data for Bayesian analysis 

The concatenated SA data usually contain more than 5 000 visi- 
bilities per channel which is not tractable for our Bayesian analy- 
sis package. The data therefore undergo a binning step using a bin 
width set to a fifth of the width of the aperture illumination function 
which leaves about 1000 visibilities per channel, evenly spaced in 
i(v-space, to be used for further analysis. The bin width was chosen 
in order to assure that enough samples populate each bin and no 
information will be lost during the process. 



4.4 McAdam - Bayesian Analysis 

Our Bayesian Analysis package, McAdam Ma rshall et al. I ( l2003h 
uses a fast sampler, muufinest Feroz et al. ( 200^^ and performs a 
joint fit of the i(v-data for the presence of a cluster imprint via a 
physical or analytical model and contaminations from undesired 
signals originating from CMB primordials, radio sources, instru- 
mental noise and source confusion. This is done in a fully Bayesian 
way; the sampling takes into account all the prior knowledge of the 
cluster model and the source environment investigated by the LA 
in order to do its simultaneous fit. 

Despite this efficient use of initial knowledge which consider- 
ably narrows down the size of the parameter space to be explored, 
the dimensionality of the problem (e.g. 40 dimensions in the case 
of ISCS 1438-1-34), sometimes requires further simplification of the 
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problem by fixing the flux and spectral index of sources whose inte- 
grated fluxes are less than 4crsA, where ctsa is the continuum noise 
level on the SA. For this analysis, we use two different models to 
describe and fit for the cluster and its imprint on the sky, as follows. 



4.4.3 Source parameters 

In both cases, the same source parameterisation is used; we fit their 
positions with a delta prior as the LA is able to measure positions 
to greater precision than the SA. All of the sources brighter than 
four times the noise level on the SA maps are fitted. 



4.4.1 /3-model 

First, w e use an isothermal fl-model ICavaliere & Fusco-Femianol 
([197^ f 1978h which has proven itself before; see e.g. 
IAMI Consortium: Hurley- Walker et al. I ( l2012h for a compari- 
son of mass derivations using an SZ yg-model and weak- l ensing 
dat a but alsolAMI Consortium: Rodrfguez-Gonzalvez et alj jlOllI 



and I AMI Consortium: Olamaie et al. I j2012l) . The model describes 
a cluster gas electron density rig which decreases with radius r 
according to 



P.(0) 



(1) 



where Pg(0) is the gas mass per electron (Pg(r) = n^ir) x 1 . 14m,„ nij, 
is the proton mass) and r^ is the core radius. We measure the total 
mass Mt.200 at radius r2oo, the radius at which the mean enclosed 
density is 200 times the critical density pciit at the cluster's redshift 
z and assuming spherical symmetry; 



47T 

Mt.200 = yr^oo(200pcrit). 



(2) 



Following lVoitldioO^ and assuming the cluster is virialised at 
this radius, we can relate a collapsing top-hat density perturbation 
model to a singular truncated isothermal sphere. This also takes 
into account the finite boundary pressure and assumes all kinetic 
energy is internal energy of the hot plasma. Hence, we calculate 
the temperature at ^200 via a mass-temperature scaling relation 



^87-^,200 = ^ f^l ' X [GMj,2mmz)]-" . 



(3) 



In summary, the data set can be fitted to a full cluster model 
with parameters 0^ = {^c,yc,fi,rc,Mj200,z,fgas.20o) and radio 
source parameters T = (x„y„So,a). A summary of the priors on 
each parameter is given in Table |4] 



4.4.2 DM-GNFW model 



Thi s parameterisation uses an analyt ical model described 
m amaie. Hobson. & Grained ( l2012h which models the 
dark mat ter halo of a cluster o f galax ies with an NFW- 
profile jNavarro, Frenk. & White I 1 1997b and a GNFW 
( iNagai. Kravtsov. & VikhlininI |2007|) pressure profile to de- 
scribe the cluster gas. This model relies on 3 sample parameters, 
which are all given at r2oo'- Mt,2oo, feas and z. We calculate 
the halo concentration parameter C200 as a function of cluster 
mass and redshift calculated using a relaxed cluster relationship 
dOlamaie. Hobson. & Graingell2012l ; lNeto et alj20d7h : 



C200 ■ 



5.26 



M- 



T,200 



lO'^/i-'Mf 



(4) 



The lower section of Table |4] shows the priors used for the DM- 
GNFW model analysis; note that the cluster's positional and red- 
shift priors are identical to thosed used in the y6-model. 



4.5 Detection of a cluster 

Using Bayesian analysis carries the great advantage that we can 
compare two hypotheses, namely the existence of a cluster (Hi) 
in a particular search area and the hypothesis of there not being a 
cluster in the same area (Hg), using the ratio R: 



R ■■ 



Pr(H,\D) ^ Pr(Z)|/fi)Pr(//i) 
Pr(Ho\D) Pr(Z)|/fo)Pr(//o)' 



(5) 



R depends only on the ratio p'^^f/J^ of the prior probabilities and 

the evidence ratio pjjpj^^j , which is an output from multinest. The 
ratio of prior probabilities can normally be set to unity but occa- 
sionally requires further consideration. Hence a formal detection 
criterion can be derived by analysing every data set twice; once al- 
lowing McAdam to fit for a cluster imprint in the observed sky and 
once prohibiting the presence of any SZ signal. Furthermore, from 
the marginalised posterior probability distributions of each param- 
eter we can extract constraints on the fitted parameters and hence 
conduct reliable parameter estimation which further enhances con- 
fidence in the detection. [Note that this R value, suitable for use 
when a cluster z is k nown, is different from the R value for b lind 
detection discussed in lAMI Consortium: Shimwell et al.l ( |2012[) .1 



5 SIMULATIONS 

In addition, to analysing the AMI SZ data, we also simu- 
late the clusters using X-ray fitted y6-model estimates from 
ICulverhouse et al.1 ( l2010h which are listed in Table |5] together with 
the simulation-derived fg^s, Mf^oo, and integrated SZ flux in AMI- 
SA channel 5. It is evident that XMJ0849+4452 in particular has a 
very low predicted fgas, which is a symptom of three problems: 

(i) Assuming isothermality with a temperature measured near 
the cluster core (excising any cooling flow), and doing the same for 
gas mass, will each i ntroduce a bias in the mass estimates at high 
radius. For example, iMroczkowskil ( l20Ilh found that the deriva- 
tion of the total mass using SZ data and the virial theorem alone 
is probably most impacted by the assumption of constant gas mass 
fraction. 

(ii) Estimating the mass internal to a radius of r25oo from X-ray 
data (using e.g. 'onion peeling') can work well. However, for in- 
struments (such as AMI and SZA) that are sensitive to larger scales, 
a problem arises along the line-of-sight of y up to a fiducial ra- 
dius r. 



(6) 



Y= f ydO. = r dl { n,{r')T,{r)Anr'^dr' . 

Ja fne& J_„ Jo 

The measured SZ flux of the telescope is related to its sensitivity 
to a certain scale on the sky rather than to a defined radius, which 
means that the actual line-of-sight contribution of Y between the 
defined radius (("2500) and ±00, is not taken out when relating Y to 
the cluster temperature and mass parameters because, for example, 
the true value of f} at large radius may be significantly different from 
the one assumed closer in. There is no easy way to correct for this. 
The SA has one advantage in these matters: because of its wider 
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Table 4. Priors for the cluster and source pai'ameters used for the Bayesian analysis of our data. We list the parameterisation of the /J-model first, then the 
source parametersation and prior distributions for each individual object which are used for both cluster models. The last section lists the alternative parameters 
and priors when using the DM-GNFW cluster model. 



Parameter 



Prior used 



Cluster Position (Xc, yc) 
Redshift (z) 
Core radius (r^ /kpc)) 
Beta (fi) 

Mass (M7.200 IMq) 
Gas fraction (/g,20o) 



Gaussian at X(;iuster5 — 60" 
(5-function on the redshift provided in Table|2] 
Uniform between 10 and 1000 
Uniform between 0.3 and 2.5 
Uniform in log space between 3.0 X lO'^ and 5 X lO" 
Gaussian prior centred on 0. 123^^son_etaL QOl^; Zhang_etaL QOlO) with tr = 0.03 



Source position fej^,) 
Source flux density (So/Jy) 
Source spectral index {a) 



A delta-function prior using the LA positions 
A Gaussian centred on the LA continuum value with a o" of 40 per cent 
A Gaussian centred on the fitted spectrum and the LA en'or as tr 



Cluster Position (Xc, jc) 
Redshift (z) 
Mass (M7.200 /Mo) 
Gas fraction (/g,20o) 



Gaussian at Xt-iuj^te^' ^ — 

60" 

i5-function on the redshift provided in Table |2] 
Uniform in log space between 3.0 X lO" and 5 X lO'^ 
Gaussian prior cenfi'ed at 0. 1 with cr = 0.02 



field of view and synthesised beam size, any measurement will give 
values of temperature and gas fraction that are averages internal to 
r2oo and hence deliver less biased (in these respects) estimates of 
the mass and temperature. 

(iii) Another issue arising from a small-radius fit to estimate 
the quantities internal to a larger radius is that any (small) er- 
ror will cause large uncertainties on quantities calculated at r2oo- 
Again, since the n^-weighted X-ray temperature is more sensitive to 
clumping and shocking than the ^-weighted SZ temperature mea- 
surement, such X-ray measurements are biased high. On the other 
hand measuring the total mass out to a radius of r^oo using this 
parameterisation would bias the temperature towards lower values 
even though we emphasise that we measure an average temperature 
on the SA. 

Ideally, one would try to measure a radial /^a, and Mtot profile 
via a joint analysis of different instruments sensitive to different 
radii. With this in mind, the following simulations and comparisons 
should be regarded as a rough estimation of cluster detectability 
rather than an accurate parameter comparison. Two targets in the 
sample (ISCS 1438-1-34 and SPJ1638-I-4039) were not simulated as 
there were no X-ray fitted values no r any SZ mass measurements 
provided in'Culverh ouse et al.l j2010l) . For SPJ1638+4039, we use 
a mass estimate at r 2oo derived from measured velocity dispersions 
jMuzzin et ai]|2009l) to give an educated guess on the detectability 
of this target. The fact that they are both clusters that have been 
successfully identified in IR but not at radio frequencies hints at a 
complex radio environment or a cluster mass that is well below the 
thermal SZ detection limit of the observing telescope. 

Our in-house package profile can simulate a patch of sky 
populated by a galaxy cluster, point sources and primordial CMB 
contribution. We use the source information from the LA raster 
observations to model the known sources in the field and add 
a population of fainter sources follow ing the IOC source counts 
lAMI Consortium: Davies et al. I lloi J), which serve to simulate 
source confusion noise, profile lets us simulate the process of ob- 
serving accurately by sampling in Mv-space according to the ar- 
ray's configuration, observing times and frequencies and adds in- 
strumental noise (assumed Gaussian) to each visibility. The simu- 
lated MV-data are then processed through the same data reduction 
and analysis pipeline as the real data. Since all of the simulations 
were generated using a /?-model, it does not make sense to analyse 
them using another model (i.e. the DM-GNFW model) as a fitted /3- 



model will always be the best match to our simulations and deliver 
the best evidence of detection. We also investigate the parameter 
constraints of each cluster to predict an estimated integrated SZ 
flux of the cluster. This analysis will contribute to estimating the 
feasibility of a real detection from the real data. 



6 RESULTS AND DISCUSSION 

For each target, we discuss the source environment, the simulation 
of the observation conducted and hence the prospects of a success- 
ful cluster detection and source subtraction. In addition, we provide 
evidence ratios and present the ID and 2D posterior marginals to 
formalise a detection. 

We present the maps of the SA observations in Fig.[T]and the 
source-subtracted maps after Bayesian Analysis in Fig.|2]for each 
target. CLJ1415H-3612 is presented first because it is the most con- 
fident detection, followed by XMJ0849-)-4452, a prime example of 
a non-detection. 



CLJ1415-H3612 

This cluster has a wl.3 mjy source located directly on top of it, a 
0.9 mJy source » 1.2' to the north and a 0.5/0.6 mJy pair » 3' to 
the north. Moreover, three other sources lie about 5' to the east 
with integrated flux densities 0.7, 0.8, 0.2 mJy respectively; the 
central source of the three could be matched in both SDSS and 
2MASS data to NGC 5529, a galaxy in a group at redshift 0.009 
jKochanek etai]|200ll) . Fig. |l(a)| shows the SA map. We note as 
well that the integrated fluxes of the three sources summed together 
are lower on the LA map than the unresolved single feature on the 
SA map in the same location. This is evidence for extended emis- 
sion, which would be plausible for a close-by galaxy or group of 
galaxies. As these sources are comparatively low in flux, we man- 
aged to successfully subtract the sources from the data using our 
fitted values and recover a 5 - cr decrement with an integrated flux 
of « 780 yuJy (Fig.[2(a)l(. 

Judging from the simulations, we do not expect to detect this 
cluster. The evidence ratio is in favour of a non-detection and the 
parameters are not constrained by the simulated data, which can 
be attributed to the fact that this simulated cluster has a small gas 
mass fraction and core radius, which leads to a small gas mass 
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Table 5. Values from /3-model X-ray fits taken from' Culverhouse et al.l j201(^) that we use for the AMI simulations, plus the AMI simulation results for the gas 
mass fraction, total mass at r2oo. and integrated SZ flux density in the 15.5-GHz channel. 



Cluster 
name 



kT 
(keV) 



Lx 



(arcsec) 



/? /ga. 



(xlO"*Mo) 



Integ. flux 

0"Jy) 



CLJ1415+3612 



RDJ09 10+5422 
XMMJ0830+5241 



6.5+ 



XMJ0849+4452 6.7^2.0 
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7 6+° '* 



io![ 
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^•'-0.2 

16+1 
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"■"■^-0.08 

0.83+oS 



10 9+" '' 
12 1+'-' 
17.9!3;0 
28.6+'o 
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0.7 0.035 
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7.281 
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Figure 1. SA maps of the cluster sample. The contours are scaled linearly in multiples of the map noise level, starting at the 3-0" level. Negative levels are 
shown with dashed lines and positive ones have solid contours. The synthesised beam is shown on the lower LHS. The positions of the sources found by the 
LA (see section l4^ in the field are displayed using crosses. 



(~ 5.266 X IO'^Mq). In complete contrast, the AMI data show a 
clear SZ detection: the model selection ratios of for each parame- 
terisation indicate a decisive detection (R > 8) which is substanti- 
ated further by the good mass constraints of each model, presented 
in Table [7] Both derived masses are substantially higher than the 
simulation predicts (Mt,200 ~ 4 x 10''* Mq), the y6-model estimates a 
mass of ^ 7.3 X 10''' Mo which agrees with « 10.4 X 10'-* Mo from 
the DM-GNFW model within 2cr. 



We show the constraints on each parameter of the cluster 
model via their ID and 2D marginalised posterior distributions for 
the /3-model (Pig.^ and (Pig.^ for the DM-GNFW model. 



XMJ0849-I-4452 

This field suffers from substantial contamination from bright 
sources, the brightest (s; 10 mjy) is located 2.75' to the south-east, 
followed by a a; 1 mJy source to the north-west about 4' away. Sev- 
eral fainter sources (w 2 mJy) also surround the pointing centre and 
are about 3'away. The simulation of this cluster shows how the cen- 
tral bright source obstructs the cluster decrement and the analysis 
of the simulation does not predict a convincing detection (R=0.05). 

From the AMI data, our evidence ratios R are < for both 
models. We present the ID posterior marginals in Fig. |4] for the 
yS-model to provide an example of what typical marginalised pos- 
terior distributions look like for a non-detection. Since the data do 
not contain any information to constrain the cluster parameters, we 
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Figure 2. Source-subtracted maps of the cluster sample observed by the SA. The contours are scaled linearly in multiples of the map noise level, starting at the 
3-0" level. Negative levels are shown with dashed lines and positive ones have solid contours. The synthesised beam is shown on the lower LHS. The positions 
of the subtracted sources are denoted by crosses. The boxes show the position of the clusters found by McAdam. 
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Figure 3. Synthesised beam of the SA for XMJ0849+4452. Contours start 
at the 6% level for both positive (solid lines) and negative (dashed lines) 
scales and increase/decrease in 3 per cent steps. 




200 400 600 800 1000 



-200 -100 





Figure 4. ID marginalised posterior distributions (black) for the parame- 
ters of XMJ0849+4452 analysed with the /3-model. The red dashed lines 
display the prior distribution for each parameter and hence show the influ- 
ence on each posterior. The green bars show the mean of each posterior 
distribution. 
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only recover the prior distributions or are heavily biased by them. 
To illustrate this point, we added the prior distributions in red in 
FiglH all of total mass, and /? are fully prior-driven in this exam- 
ple. 



ISCS1438H-34 

This observation is strongly contaminated by two bright sources 
(for a map, see Fig. |l(c)^ , one with an integrated flux of about 
38 mjy, approximately 8.5' to the east (14 38 49,+34 16 00) and the 
another one, with integrated flux density of about 18mJy, a 9' to 
the south (14 38 11, -1-34 05 06). Both sources are well embedded 
within a crowded environment and produce significant sidelobes. 
Although none of them are classified as extended on the LA, the 
parameter space to be explored in order to model them accurately 
is too large for McAdam. If we fix all the source fluxes and spectral 
indices using delta priors from the LA measurements, apart from 
a few sources closest to the phase centre and deemed to cause the 
highest nuisance to the data, we reduce the dimensionality of the 
Bayesian fit, which improves convergence on a best fit in the anal- 
ysis. However, we are sill not able to recover an SZ decrement nor 
any parameter estimates from the data. The model selection crite- 
rion strongly favours a non-detection. 

RDJ0910-I-5422 

Two sources with integrated fluxes of 36 and 17 mJy respectively 
on the LA map produce important sidelobes which swamp the en- 
tire map (Fig. |l(d)y The two sources form one extended feature on 
the SA map 7' to the north-west of the phase centre. 

We do not expect this cluster to be detectable; Table |5] indi- 
cates a relatively low-mass cluster with a low integrated flux esti- 
mate. In the simulation however, we do manage to successfully re- 
cover some SZ flux despite a moderately successful source fit and 
subtraction. Judging from the map of the simulated data, shown in 
Figure[5] We expect a negative peak flux from the cluster of « -400 
yuJy in the real data, which is less than 1 per cent of the brightest fea- 
ture in the map and will be a potential issue for the dynamic range 
of the telescope due to correlator errors and residuals left from the 
source subtraction. This makes a successful SZ signal recovery and 
source-subtraction from real data very doubtful and indeed the de- 
rived detection corroborate the hypothesis of a non-detection. 



SPJ1638-(-4039 

This target has a relatively clean source environment (Fig. |l(e)[ (; 
the closest source from the phase centre is » 5' away and has an 
integrated flux of about 1 mJy. If we use the estimated mass at 
^200 derived from m easured velocity dispersions of SPJ1638+4039 
jMuzzin et al.l2009l) . Mom = (2.4 ± 1.8) x 10"* Mo and the fact that 
this cluster is at high redshift (z = 1.2), we do not believe that this 
cluster is able to produce an SZ signal which is above our detection 
threshold. A simulation of the cluster's imprint on an empty sky, 
using the DM-GNFW model and assuming a gas mass fraction of 
0.123 gives an integrated flux estimation of w fijy , which would be 
a 2-cr detection. In the SA data, the source subtracted map (Fig. 
|2(e)^ does not show any evidence of SZ flux and reflected by a 
negative model selection parameter for both models in Table [6] 




09 12 00 11 30 00 10 30 00 09 30 

Right ascension (J2000) 

Figure 5. Map of the simulated obsevation of RDJ09 10+5422. The con- 
tours are scaled linearly in multiples of the noise level, starting at the 3(T- 
level. Dashed lines depict negative flux levels and solid contours show pos- 
itive levels. 



Table 6. Evidence ratios R for each cluster and different parameterisations. 
Also provided are the evidence ratios for the simulations. 

Cluster /3-model DM-GNFW model simulation 



CU1415+3612 8.6 8.0 -0.1 

XMJ0849+4452 -0.64 -0.2 0.05 

ISCS1438+34 -6.1 -5.4 no data 

RDJ0910+5422 -3.2 -2.9 0.6 

SPJ1638+4039 -1.1 -1.9 no data 

XMJ0830+5241 3.0 3.2 4.3 



XMJ0830-I-5241 

Fig. [1(0] shows the SA map of XMJ0830+5241, which has a 
moderately crowded source environment. However all the closest 
sources, within a distance of 1 synthesised beam, are relatively low 
in flux (1 mJy and less). There is another source about 3.5' to the 
south-west, which has an integrated flux density of » 2 mJy. Fur- 
thermore significant decrements can be seen on the map, prior to 
any source subtraction which cannot be assigned to the sidelobes 
of any radio sources in the field (see Fig. [3] for a reference syn- 
thesised beam). Simulating this observation with the LA measured 
source fluxes and spectral indices as well as the cluster parame- 
ters from Table [5] strongly favours a detection (R > 30), which is 
already clear by mapping the simulated uv-data in which the clus- 
ter's imprint appears with more than 6 times the noise level of the 
observation. Analysing the SA data, the evidence ratio for a detec- 
tion is substantial for both parameterisations (R >3) although less 
strong than in the case of CLJ1415+3612. After source subtrac- 
tion, the integrated flux of the the cluster's imprint is a 600 /jJy . 
The posterior marginals, shown in Fig. [7] for the /3-model and in 
Fig. |9} for the DM-GNFW model, show good constraints on the 
mass for both parameterisations; they agree with each other within 
one cr and indicate a relatively low-mass cluster of » 4 x 10'"* M©. 
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Figure 6. ID and 2D marginalised posterior distributions for 
CLJ1415+3612 (/?-model). Tlie green crosses and lines show the mean 
of the respective marginal distribution. The contours on the 2D marginals 
indicate the areas enclosing 68% and 95% of the probability distribution. 
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Figure 8. ID and 2D marginalised posterior distributions for 
CLJ1415+3612 (DM-GNFW model).The green crosses and lines show 
the mean of the respective marginal distribution. The contours on the 2D 
marginals indicate the areas enclosing 68% and 95% of the probability 
distribution. 




-50 50 -100 100 500 1000 0.5 1.5 2.5 2 4 6 8 1 2 



Figure 7. ID and 2D marginalised posterior distributions for 
XMJ0830+5241 (yS-model). The green crosses and lines show the 
mean of the respective marginal distribution. The contours on the 2D 
marginals indicate the areas enclosing 68% and 95% of the probability 
distribution. 
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Figure 9. ID and 2D marginalised posterior distributions for 
XMJ0830+5241 (DM-GNFW model). The green crosses and lines 
show the mean of the respective marginal distribution. The contours on the 
2D marginals indicate the areas enclosing 68% and 95% of the probability 
distribution. 



Table 7. Mean values and their 68% confidence Hmits of Mt.200 for each 
cluster and parameterisation. 

Cluster Mt.20o/Mo Mj2oo/Mq 

name CS-model, SA data) (DM-GNFW, SA data) 

CLJ1415+3612 7.27+] X 10''' 10.40^^,, x lO'" 

XMJ0830+5241 3.56+] ['J x 10''' 4.66+] :J| x lO'" 



7 LA SZ MEASUREMENTS 

For our two SA SZ detections, XMJ0830+5241 and 
CLJ1415+3612, we do additional follow-up observations 
with the LA. As both clusters are at high redshift (0.99 and 1.03 
respectively) we aim to investigate whether these clusters' angular 
extent would be better geared to the baseline range of the LA 
and hence improve parameter estimations. Note that our Bayesian 
analysis expects the cluster's imprint to be an extended feature 
and hence have a different profile as a function of j(v-distance 
than a point source. The trade-off on the LA is the more sparsely 
filled Mv-coverage which will diminish the sensitivity to the SZ 
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effect. From our SA analysis, we find r2oo ~ 2.3'and r2oo ~ 2.9 'for 
CLJ1415+3612 and XMJ0830+5241, respectively. As the angular 
resolution of the LA is » 30" using all of its baselines, these two 
clusters will be resolved by the LA. 

We note that the maps agree with what we expect; the 
LA map of XMJ0830+5241 (Fig. [70(a)l ( show a 3<t decrement 
which agrees with the SA position. After source-subtraction (Fig. 
|10(b)^ , we manage to recover a 5cr decrement. The deep follow-up 
CLJ1415-I-3612 is swamped by the « 1.3mJy located on top of the 
cluster position, as shown in Fig. |10(c)| After source-subtraction, 
there is a small, just over 3cr, decrement in the source-subtracted 
map whose position agrees with the SA observation. 

We cannot do any Bayesian analysis of the LA data with our 
current models, which requires application of the virial theorem 
internal to raoo- 



8 CONCLUSIONS 

AMI has observed a d eclination-limited subset of six of the 
ICulverhouse et al.1 ^201(1) sample. 

(i) Three of these targets (ISCS 1438+34, RDJ0910+5422, 
XMJ0849-I-4452) suffered from heavy radio source contamina- 
tion which inhibited a detection, one other (SPJ1638+4039) was 
deemed to produce a thermal SZ effect which is too small to be 
detected. 

(ii) Two were confirmed detections (CLJ1415+3612, 
XMJ0830-I-5241) with the SA, with Bayesian evidence ratio 

of a 8 and » 3 respectively. 

(iii) In the subset of six, AMI detec ts in SZ the same c lusters 
(CLJ1415+3612, XMJ0830+5241) as ICulverhouse et al.l i2010h 
with the SZA. The AMI data and analysis, however, return pa- 
rameters internal to rim (~ virial radius). For CLJ1415+3612 and 
XMJ0830-I-5241 respectively, our yS-model, route finds Mt,2oo = 
7.3 ± 1.8 X 10''*Mo and Mt.2oo = 3.6 ± -I- 1.1 x 10''* Mq, while 
our DM-GNFW model finds Mt,2oo = 4.7 ± 1.4 x IO^Mq and 
Mt,2oo = 10.4 ± 2.5 X 10"* Mo. 

(iv) In using the X -ray data internal to ^2500 from 
ICulverhouse et all j2010l) to simulate what AMI should see 
at larger radius, we have highlighted three likely biases that an SZ 
measurement will produce with a low-radius X-ray measurement. 

(v) We find higher-resolution SZ effects in the follow-up obser- 
vation with AMI LA of CLI1415+3612 and XMI0830+5241. 
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Figure 10. Maps of the deep single pointed observations of XMJ0830+5241 (top) and CLJ1415+3612 (bottom) on the LA before (left) and after (right) source 
subtraction. The mapping was done using a wv-taper of 1.2 ki. Crosses show the positions of subtracted sources and the box indictates the position of the 
cluster found by McAdam (see text for further explanations). Contours scale in multiples of the noise level, starting at 3 — o"; negative contours are displayed 
with dashed lines while solid ones denote positive ones. The synthesised beam is displayed as an inlaid on the lower left-hand side of the map. 
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